Scaling theory of two-dimensional metal-insulator transitions 
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We discuss the recently discovered two-dimensional metal-insulator transition in zero magnetic field 
in the light of the scaling theory of localization. We demonstrate that the observed symmetry 
relating conductivity and resistivity follows directly from the quantum critical behavior associated 
with such a transition. In addition, we show that very general scaling considerations imply that any 
disordered two dimensional metal is a perfect metal, but most likely not a Fermi liquid. 
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In an influential 1979 paper [Q, simple scaling argu- 
ments were put forward which have shaped much of our 
thinking about metal-insulator transitions (MITs) in dis- 
ordered systems. Focusing on noninteracting electrons, 
the authors demonstrated that in two dimensions (2D) 
even weak disorder is sufficient to localize the electrons, 
and no true metallic behavior is possible at T = 0. Ever 
since, it has been widely believed that two is the lower 
critical dimension for MITs in disordered systems and 
this prompted much theoretical activity based on ex- 
pansions around dimension d — 2. The scaling ideas 
were quickly extended to interacting systems, but despite 
many years of painstaking effort, no satisfying theoretical 
picture has emerged for the interacting case. Neverthe- 
less, the belief that all the states are localized in d = 2 
has remained largely unquestioned. 

Recently, a remarkable experiment performed on a 
2D electron gas in zero magnetic field suggested that the 
conventional picture might be incomplete. In this work, 
fairly convincing evidence was presented which supports 
the existence of a true MIT in 2D. Despite the beauty 
and elegance of the data, this work has met considerable 
resistance, largely due to its apparent conflict with the 
scaling theory of localization as well as its uniqueness. 

The major assumption of was based on an earlier 
idea of Thouless [|j about the length (scale) dependence 
of the conductance. In [|l| it is asserted that the "beta 
function" P{g) = d[log((7)]/(i[log(L)] is a function of the 
conductance g itself, but not an explicit function of the 
length scale L. Now, P{g) is known in the two limits of 
very large and very small disorder, and it is reasonable to 
assume that it is continuous (smooth) in between. From 
Ohm's law, P{g) = {d — 2) for g -\-oo, while for g 
small, one expects exponentially localized states, so that 
f3{g) — *■ — oo in this limit. Since the metallic behavior 
is possible only for /3(g) > 0, the form of P{g) at large 
g is sufficient to determine the stability of the metallic 
phase. In particular, for noninteracting electrons, this 
has a form P{g) = {d — 2) — l/g -|- ... |^, indicating that 



the metallic phase is unstable in d < 2. 

However, it is important to emphasize that the last 
step in this analysis is valid only for noninteracting elec- 
trons; the scaling theory is only complete in the absence 
of electron-electron interaction. In the following, we show 
that the existence of a 2D transition does not contradict 
any general scaling principles. However, we argue that 
within the scaling theory such a transition has a number 
of unusual features, many of which are already apparent 
in the existing data. 

No fundamental principle requires that f3{g) be mono- 
tonic, or in fact negative in 2D when g is large. Instead, 
if we assume that the leading behavior at large g is 

/3(g) = (d-2) + A/g" + ..., A>0, (1) 

then it is clear that the beta function for all d > 2 has to 
change sign at some finite g = gc, leading, in particular, 
to a metal- insulator transition in 2D. Alternatively, the 
existence of a MIT in 2D requires that /3(g) > at 5 
large, where the exponent a simply parameterizes how 
the system approaches the metallic limit as 5 — > 00. 

This idea is not new. Such a scenario is realized for 
non-interacting electrons in the presence of spin-orbit 
scattering In this case, a — 1 and a MIT is expected 
in 2D. However, this behavior has never been taken too 
seriously in the context of real 2D systems because it is 
known that adding interactions to the spin-orbit univer- 
sality class results in reversing the sign of the first quan- 
tum correction and all the states again end up local- 
ized in 2D. A more complicated situation was found in 
the absence of spin-orbit and time-reversal breaking per- 
turbations. There, an interaction driven enhancement 
of the conductivity was found at weak disorder which 
was expected to overwhelm the localizing effects of coher- 
ent backscattering (weak localization) However, the 
analysis revealed that the effective interaction strength 
diverges upon scaling, making it difficult to determine 
what will actually happen at long scales or low temper- 
atures m. Nevertheless, from what is known it appears 
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that the relevant quantum corrections due to interactions 
can only enhance the conductivity at low temperatures, 
so that the possibility of a 2D metal remains open. 

At present, the abovementioned problems in the in- 
teraction case make it difficult to determine the form of 
the leading quantum corrections even at weak disorder. 
Rather than elaborating on this, we shall assume that 
interactions lead to a positive (3 at g large, so that the 
2D metal exists. The first question that one can address 
is what is the low temperature form of the conductivity 
in the metallic phase? Within the scaling theory this is 
determined only by the form of f3{g) at g large. Assum- 
ing a form of Eq. (1), and integrating the beta function 
from £ (mean free path) to L (sample size), we find 



5(L) = [.go" + ^«log(L/^)]i/". 



(2) 



If we consider an infinite sample at finite temperature, 
then Leff ^ T^^^^. Here z is the dynamical exponent 
associated with the metallic {g = oo) fixed point. One 
expects z = 2 in a conventional diffusive regime, in the 
absence of dangerously irrelevant variables Q. We thus 
find that at low temperatures in the metallic phase 



g{T) ^ logi/"(To/T). 
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In other words, the conductance will diverge^ i.e. the 
resistance will vanish at T = in the metallic phase, al- 
beit in a weak, logarithmic fashion. This should be true 
at sufficiently low temperatures throughout the metallic 
phase, i.e. the temperature dependence should have a 
universal form, as in Fermi liquid theory. However, we 
emphasize that this is not a Fermi liquid, not only be- 
cause of its unusual temperature dependence, but more 
fundamentally, since we expect the nature of these elec- 
tronic states to change completely if the interactions are 
turned off: they would localize. 

Another experimentally relevant question is what is the 
temperature dependence of the conductance in the quan- 
tum critical region associated with the metal-insulator 
transition. Within scaling theory, this question can be 
answered very precisely, as follows. In our formulation, 
the transition occurs a,t g = gc where the beta func- 
tion changes sign, i.e. (}{gc) = 0. We recall that deep 
in the insulating phase, the beta function is logarithmic 
in g |jl[. It is thus plausible to introduce the quantity 
t — log(g/gc) as a natural scaling variable. To deter- 
mine the critical behavior it is sufficient to consider the 
linear approximation to the beta function near its zero 
at t = 0. The slope of l3{t) at t — determines the 
correlation length exponent v. 
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t + 0(t2) 



rf(logL) 

By integrating this equation from £ to L, we find 
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where to = log(go /gc) is determined by the starting value 
go of the conductance at scale £. In the critical region, 
where we start with a go very close to g^ to ~ {go — 
9c) 1 9c ^ ^iT-j where n is the carrier concentration, and 
Sn = (n — nc)/nc measures the distance to the critical 
point. Then from Eq. (5), we have for the conductance 



g{L) — gcGxp ASn{L/£) 



(6) 



At non-zero temperature, the length scale L is deter- 
mined by the temperature through the dynamical expo- 
nent as T and we obtain the temperature depen- 
dence of the conductance 

g{6n,T) = g.exp {sgn{Sn)A[To{Sn)/T]^) , (7) 

where we have defined a crossover temperature To corre- 
sponding to the inverse correlation time as 



ToiSn) ^ \Sn\ 



(8) 



Here, A is an unknown dimensionless constant of order 
one. Let us define the scaled conductance as g*{5n,T) — 
g{Sn,T)/gc. From Eq. (7), we then immediately find a 
striking symmetry relating the conductance on the metal- 
lic side {Sn > 0) to the resistance on the insulating side 
{6n < 0) of the transition at g* ~ 1 (i.e. at Sn = 0): 



g*i6n,T)^l/g*{-Sn,T). 



(9) 



Remarkably, precisely this behavior is clearly seen in the 
experiments |2j, as emphasized in Rcf. 

In related work pO[ | , a similar symmetry for the depen- 
dence on electric field was found for the same 2D system. 
If one uses the conventional electric field scaling hypoth- 
esis 11 1, the preceding arguments carry over with the 
electric field E replacing T and 1 + z replacing z, so that 
this symmetry too appears as a consequence of the scal- 
ing argument we have presented. In fact, the behavior 
holds whatever the exponent relating the electric field to 
a characteristic length. 

In order to fully appreciate the significance of these 
findings, we should carefully qualify them and comment 
on the range of their validity. 

(1) In the authors point out a similarity of this 
symmetry with the one observed in the context of the 
quantum Hall liquid-to-insulator (QHI) transitions [[l2[ . 
However, we emphasize that the theory described in |12| ] 
is fundamentally different from that of the present pa- 
per. The symmetry found in a 2D MIT is restricted 
to the quantum critical region, which is defined by the 
crossover temperature scale, i.e. it is expected to hold 
only for T > To{Sn). In fact, the experiments completely 
confirm this expectation, as the authors themselves point 
out 1^. A careful examination of the data in ||^ reveals 
that the temperature at which a departure from sym- 
metry is observed is T « To{Sn), with Tg shown in the 
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inset of Fig. 3 of The same feature is clearly seen 
also in the earlier data of where in Fig. 4, the resis- 
tance is plotted as a function of the reduced temperature 
T = T/T()(Sn). Here, the obvious symmetry of the re- 
sistance and the conductance on the respective metallic 
and insulating sides of the MIT is seen and it is violated 
below T « To{Sn), as expected. The duality observed in 
QHI is thought to be a charge-flux duality that holds for 
the full Chern-Simons Lagrangian, ignoring the roles of 
irrelevant operators, cutoffs, or disorder. If this is so, it 
has little resemblance to the symmetry in MIT. On the 
other hand, it may be that the data in QHI transition 
can be understood from the present perspective, where 
{B — Be) oc [l/i/ — l/vc) is the control parameter that 
tunes the system through a quantum critical point; here 
B is the magnetic field, and ly is the filling fraction. 

(2) We emphasize that the above results were derived 
by linearizing the beta function around the transition 
point g = gc- The symmetry is then simply a conse- 
quence of the smoothness of the beta function in the crit- 
ical region. However, without further knowledge about 
the form of the beta function, the result is strictly speak- 
ing valid only for Sn <^ 1, i.e. it only defines the form 
of the leading high-temperature correction in the quan- 
tum critical regime. In other words, to be consistent, one 
should keep only the leading term of the expansion of the 
exponential of Eq. (7): 



gi6n,T) Kgc\l + ASn/T^^/"^^ 
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This expression is generally expected to be valid only for 
\g{T) — gd/gc ^ 1- On the other hand, if we assume that 
the linearized expression for the beta function [Eq. (4)] is 
a good approximation over an appreciable conductance 
range gmin < gc < gmax , then the full exponential behav- 
ior is valid in that range. The experimental data reveal 
that the symmetry is found in a much broader range than 
allowed by the form of Eq. (10). In fact, the data clearly 
display an exponential temperature dependence as in Eq. 
(7), which is particularly striking in the metallic phase, 
where the resistance is found to drop almost by an order 
of magnitude within the "symmetry" regime. The de- 
viations from leading behavior are more clearly seen in 
Fig. (2b) of |9j, where the conductance and as well as the 
resistance are plotted at a fixed temperature, for a range 
of carrier concentrations. The two curves show perfect 
mirror symmetry which holds over an extended range of 
concentrations, where the conductance shows a dramatic 
exponential dependence on concentration, in agreement 
with Eq. (7). This should be contrasted with the leading 
correction given by Eq. (10), where the density depen- 
dence is linear, as expected only for Sn I. 

(3) We conclude that the experimental results of 
P,p^p|] go beyond just confirming the expectations based 
on general scaling arguments in the quantum critical re- 
gion. They also provide striking evidence about the form 



of the beta function in the critical region. In particu- 
lar, they indicate that in a wide range of conductances 
1/4 < g/gc < 4 the beta function is well-approximated 
by the linear expression in t = ^og{g/gc) as in Eq. (4). 
How can we rationalize this finding? Deep in the insulat- 
ing regime, {g <C gc) the beta function is exactly given by 
(3{g) ~ log(g). The above experimental result can thus 
be interpreted as evidence that the same slow logarith- 
mic form of the beta function persists beyond the insu- 
lating limit well into the critical regime. Further support 
for this idea can be found in the experimental findings of 
Hsu and Valles Q on ultrathin films. While no MIT was 
found, their data yields a beta function which is close to 
linear in log{g) ("strong localization") all the way to the 
standard crossover conductance g ~ e^/2/i. If it is gener- 
ally true, this feature could be used as a basis of approx- 
imate calculations of the critical exponents. In contrast 
to the well-known 2 + e expansion, here one would try to 
obtain the form of the beta function in the critical region 
by an expansion around the strong disorder limit. 

(4) There are physical reasons which support our con- 
jecture that the beta function must be a logarithmic func- 
tion close to gc, if gc happens to fall in the regime of 
strong disorder, which the experiments discussed here 
indicate. The assumption of the scaling theory is that 
the conductance g is the only important parameter at 
large distances and low energies. That is, the distribu- 
tion of conductances of a hypercubc of linear dimension 
L is so sharply peaked that the most probable value is 
the only value, hence the mean value. This assumption 
can be justified in the limit of weak disorder, where the 
sample to sample fluctuations of the physical properties 
are weak, but not when the disorder is strong and the 
distribution of the conductances is broad. At strong dis- 
order, it is more plausible to assume, if we still want 
to describe the problem with a single scaling variable, 
that the distribution is sharply peaked when described 
in terms of \ogg [p^ , or more generally in terms of a 
function log[0((7)], where 4'{g) is a smooth function of g 
that tends to b-i/g -I- 6o -I- big + . . ., as g ^ 0, but tends 
to oq + ai/g + a^lg^ ■ ■ ■, &s g —> oo. At strong disorder, 
it is natural to believe that the Green's function between 
two points can be computed in a hopping parameter ex- 
pansion |ll5|| . It is plausible that the impurity averaged 
square of the Green's function for p hops on an effective 
lattice of spacing a larger than the mean free path is 
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If g > gc, the perturbation theory does not converge as 
p — !■ oo; this signals a failure of the localization assump- 
tion. When g < gc, this is not so, and we can define a 
localization length. However, for a finite system of linear 
dimension L, one can meaningfully use this equation for 
all g, because p is restricted to be less than {L/a). This 
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allows us to derive the beta function, as follows. From 
Eq. (11), the physical localization length is 



2a 



log [<P{9)l<P{9c)] ' 



(12) 



If we keep the physical f fixed as we vary the lattice 
spacing, that is, set d£_/da — 0, we get 
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<t>{gc 



(13) 



To check that this is correct, consider the limit gc — 0. 
We find that 



l^^ (1 + 1—9 

a log a 0-1 



•)log 



(14) 



which is precisely the locator expansion ||l|,y_5[. The dif- 
ficult task is to be able to compute 4>{9) from a micro- 
scopic model. Note, however, that local electron-electron 
interactions should not pose any special problems in the 
locator expansion. 

(5) The symmetry associated with the quantum criti- 
cal region is specific to 2D systems. In particular, stan- 
dard arguments suggest that in general dimension, the 
conductivity as a function of reduced concentration and 
temperature should assume the scaling form 



a{5n,T) - t'^'^~'^^'' JiSn/T^/"'). 



(15) 



Here, f{x) is a universal scaling function such that 
/(O) — 1, in the absence of dangerously irrelevant vari- 
ables. As a result of the extra temperature prefactor 
present when d > 2, the conductivity will vanish at the 
transition {5n = 0, T = 0), ruining the refiection symme- 
try of the quantum critical region. This fact would make 
it very difficult, if not impossible, to extract the infor- 
mation about the details of the beta function from the 
temperature dependence in the quantum critical region 
when d > 2. 

(6) The samples showing the 2D MIT are distinguished 
from those in previous work, as pointed out by the au- 
thors of 1^, by the fact that the density n is so low 
that the Coulomb interaction U oc \fn is immense, al- 
most an order of magnitude greater than the Fermi en- 
ergy Ep (X n. Furthermore, for small n, high mo- 
bility fi is required to reach the critical conductance 
9c = rice^ ~ e^/3/i. This also leads to the conclusion 
that kpl < 1 for the samples showing MIT. Thus it is 
misleading to assume that because the mobilities are high 
they are weakly disordered . These conditions both lead 
to the conclusion that the one-electron picture of ^ can- 
not be expected to be valid and may explain why the 
MIT is not observed in some samples. 

In summary, we have discussed the recently discov- 
ered 2D metal-insulator transition in zero magnetic field 
in the light of the scaling scenario for localization. We 



have shown that very general scaling considerations im- 
ply that any disordered 2D metal is a perfect metal, but 
most likely not a Fermi liquid. In addition, we have 
demonstrated that the observed symmetry relating con- 
ductivity and resistivity follows directly from the quan- 
tum critical behavior associated with such a transition. 
Furthermore, the fact that this symmetry is found over 
an extended range of conductances provides important 
information about the form of the beta function in the 
critical region. 
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